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In ferromagnets, the spontaneous magnetization bears the Hall effect through the relativistic spin- 
orbit interaction. Similar effects also occur in thermoelectric and thermal transport phenomena. 
Their mechanism, if it is of the intrinsic or extrinsic origin, has been controversial for many decades. 
We present a unified theory of these Ifall transport phenomena in ferromagnetic metals with dilute 
impurities at the zero temperature, in terms of a fully quantum-mechanical transport theory for 
multi-band systems with the self-consistent T-matrix approximation. This theory becomes exact 
with a single impurity and is appropriate for treating the dilute limit of the impurity concentration 
?^imp• With the Fermi energy Bp and the spin-orbit interaction energy E^o being fixed {Ep > -Eso), 
three regimes and the associated two crossovers are found in the anomalous Hall conductivity Oxy 
as a function of riimp that controls the longitudinal conductivity g^x ■ (i) In the superclean case with 
the relaxation rate h/r < Ui^npEsoD, the skew scattering arising from the vertex correction yields a 
dominant contribution that is inversely proportional to riimp, where Uimp is the impurity potential 
strength and D is the density of states. With increasing fi/r, this extrinsic skew-scattering contri- 
bution rapidly decays, (ii) In the moderately dirty regime UimpEsoD < fi/r < Ep, Oxy becomes 
insensitive to the scattering strength because of the intrinsic dissipationless topological Berry-phase 
contribution. It is resonantly enhanced to the order of the quantization unit of conductance, when 
an accidental degeneracy of band dispersions around the Fermi level is lifted by the spin-orbit in- 
teraction. Further increasing ft/r, another crossover occurs to (iii) the scaling regime of Oxy oc Uxx 
with </5 ~ 1.6, which has recently been verified by experiments on a wide class of ferromagnets. 
Similar behaviors also appear in the temperature-linear coefficient of the themal Hall conductivity 
hixy. The thermoelectric Hall conductivity axy strongly diverges in the clean limit when the Fermi 
level crosses edges of the avoided-crossing, which may be observed by careful experiments. With 
increasing h/r, there occurs an interference beteween positive and negative contributions to axy, 
which often leads to a sign change and obscures similar crossovers in the anomalous Nernst effect. 

PACS numbers: 72.15.Eb, 72.15.Lh, 72.20.My, 75.47.-m 



I. INTRODUCTION 

The Hall effect is a fundamental transport phe- 
nomenon in solids that an applied electric current in- 
duces a transverse voltage drop, or an applied electric 
field produces a transverse current^ii^i^. In conventional 
semiconductors and metals, this Hall current linear in a 
weak magnetic field H offers a mean to probe an effec- 
tive carrier number through the normal Hall coefficient 
Rh- In contrast to this normal Hall effect driven by the 
Lorentz force, a spontaneous magnetization as well bears 
the Hall effect in ferromagnets^., where the relativistic 
spin-orbit interaction is indispensable for connecting the 
spin polarization with the orbital motion of electrons. 
This phenomenon, i.e., spontaneous or anomalous Hall 
effect^ ''^ , has been one of the most fundamental and in- 
triguing issues in condensed-matter physics. 

Early experimental works led an empirical relation of 
the Hall resistivity p^y to the weak applied magnetic field 
and the spontaneous magnetization Af ^ both along 
the z direction; 

P:,y « RhH' + 4^i?sM^ (1) 
with Rs being called the anomalous Hall coefficient*^. 



Similar spontaneous or anomalous effects are also 
found in thermoelectric and thermal transport phe- 
nomena as the anomalous Nernst-Ettingshausen effect 
and the anomalous Luduc-righi effect, respectively. In 
spite of the intensive and extensive studies for many 
decades^i^iii^'^iiSiiiiii^ii^, a long standing debate on the 
mechanism has not been resolved yet. The keen issues 
are roles of scattering and the associated relaxation and 
dissipation. In particular, the intrinsic vs. extrinsic 
mechanisms and the associated scaling behaviorsi^ of 
the anomalous Hall effect have attracted revived inter- 
est because of the fundamental importance of the dissi- 
pationless and topological nature of the intrinsic mech- 
anism, which penetrates the whole debates on this is- 
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In a recent Letter—, we have presented a unified theory 
of the anomalous Hall effect, fully taking account of both 
the intrinsic and the extrinsic contributions on an equal 
footing. Now, main aims of this paper are (i) to pro- 
vide a comprehensive description of the unified theory 
of the anomalous electric, thermoelectric, and thermal 
Hall transport coefficients with some important details 
of the formalism and the calculation procedure, (ii) to 
explain key experimental observations in various ferro- 



2 



magnetic metals, including the magnitudes of axx and 
and their scaling relations, by means of the classi- 
fication into three regimes revealed by the theory, and 
accordingly (iii) to resolve the long standing controversy 
on the mechanism. 

The dissipationless and topological nature involved in 
the Hall effect has been highlighted by the discovery of 
quantum Hall effect^^ in two-dimensional (d = 2) disor- 
dered electron systems under a strong magnetic field. For 
the Bloch electrons in perfect crystal, the Hall conductiv- 
ity is expressed by the Thouless-Kohmoto-Nightingale- 
Nijs (TKNN) formula3^, 

a™N ^ _,^^.,e2ft^ j _|^5l(p)/(s„(p)) (2) 

with the electronic charge — e (e > 0), the Planck con- 
stant h = 2T:h, the Fermi distribution function /(e), and 
the anti-symmetric tensor eij£. We have introduced the 
eigenenergy £n{,p), the Berry-phase connection 

a-nip) ^ i{n,p\'Vp\n,p), (3) 

and the Berry-phase curvature 

bn{p) = Vp X a„{p) (4) 

of the generalized Bloch wave function \n,p) with the 
band index n and the Bloch momentum p. Each band 
is characterized by a topological integer called the Chern 
number 

The sum of C„ over the occupied bands determines the 
integer v (Chern number) for the quantization of the Hall 
conductivity axy — ve^ /h. Then, in ideal cases when the 
Fermi level is located within an energy gap, the longitudi- 
nal conductivity Gxx vanishes and the Hall conductivity 
axy is quantized in a unit of e^/ft. = 3.87 x 10^^ 57^^. This 
Berry-phase effect has been incorporated into the adia- 
batic semi-classical wave-packet equations for the Boltz- 
mann transport theory^^. 

Historically, the dissipationless thermodynamic Hall 
current was first discussed by Karplus-Luttingcr"'. They 
initiated an intrinsic mechanism of the anomalous Hall 
effect in a band model for ferromagnetic metals with the 
spin-orbit interaction. Recognizing that the inter-band 
matrix element of the current operator plays a key role, 
they derived a generic expression for the band-intrinsic 
contribution to the anomalous Hall conductivity, which 
is independent of the scattering rate. This accounts for 
the experimentally observed scaling relation for the re- 
sistivity tensor i.e., pxy oc p^^ or equivalently, Uxy 
being constant. They also performed a perturbation ex- 
pansion in M and the spin-orbit coupling ^ to derive the 
empirical law given by Eq. ([T]). However, there are two 
drawbacks in the theory. 



Firstly, the perturbation theory for Uxy in terms of 
the spin-orbit interaction energy i?so which is 

usually even small compared with the bandwidth or the 
Fermi energy Ep except in the / electrons, can not cap- 
ture a topological nature involved in the intrinsic anoma- 
lous Hall effect. Recently, it has been recognized that 
the Karplus-Luttinger's general expression for the band- 
intrinsic contribution actually coincides with the TKNN 
formula given by Eq. Q ^'^i^^ , and that each band con- 
tains a finite Chern number—, as in the quantum Hall 
systems. Without the spin-orbit interaction, the Hamil- 
tonians describing the majority and minority spin bands 
are decoupled in the band theory. Then, the fact that 
the Hamiltonians and the Bloch wavefunctions are real 
requires the accidental degeneracy of band dispersions in 
the three-dimensional Brillouin zone'^^. Turning on the 
spin-orbit interaction, this condition no longer holds and 
then the accidental band crossings are avoided, leaving 
a small energy separation of the order of -Bso- Namely, 
the spin-orbit interaction plays a crucial role in avoiding 
a crossing of band dispersions at a certain momentum po 
(see Fig. [T]). This avoided-crossing of band dispersions is 
accompanied by a transfer of Chern numbers among the 
two-dimensional bands. This phenomenon called "parity 
anomaly" in (2 + 1) dimensions has a non-perturbative 
nature^: axy exhibits a discontinuous jump by e^//i as 
Eso continuously changes its sign. This points to an 
importance of the avoided-crossing of band dispersions 
near the chemical potential. Therefore, the nontrivial 
topological structure in the Bloch wave functions of ferro- 
magnets is not captured by the perturbative treatment^ 
of the spin-orbit coupling ^ leading to the empirical law 
given by Eq. ^ with Rg oc 

This picture based on the "parity anomally" has been 
supported by recent first-principles calculations. When 
the Fermi level is located around such an avoided-crossing 
of dispersions, as found in recent ab inito calculations for 
SrRuOj2i., the bcc Fe^^, CuCr2Se4-xBr^, Co^, and 
Ni^, the magnitude of a'^™^ is resonantly enhanced; 
Cxy ^ /ha ~ 10'^ f2~^ cm~^ with the lattice constant 
a « 4 A3i^, which can be regarded as an nearly quan- 
tized axy in each two-dimensional momentum plane or pz ■ 
This resonant enhancement of axy without any small fac- 
tor of Eso ~ S,M means that the perturbation expansion 
in fails when the Fermi level is located within the en- 
ergy range of the avoided-crossing of band dispersions. In 
the metallic system, there appear many avoided-crossings 
and/or more complex structures near the Fermi level, 
which may lead to a complex behavior of axy as a func- 
tion of the chemical potential, the magnetization, and 
the crystal structure analogous to the quantum chaos, as 
actually found in first-principles calculations^-^ '^^i^^ . In- 
terference among the contributions from different bands 
and/or different momentum regions may often reduces 
the magnitude of axy, but its variation is of the order of 

/ha in both calculations and experiments on SrRuO^^. 

Secondly, scattering events extrinsic to the band struc- 
ture were completely ignored in the Karplus-Luttinger's 
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FIG. 1: The simplest examples of (a) an accidental crossing of 
two band dispersions at a momentum Pq and (b) an avoided- 
crossing with a splitting of the dispersions by 2Eso at this 
momentum region. 



theory''. In fact, adiabatic semi-classical Boltzniann 
transport analyses^i^iiSiiiiii^iii^, which have also been 
taken over to the extrinsic scenario for the spin Hall 
effecli^, have revealed that the spin-orbit interaction 
in the impurity potential produces the anomalous Hall 
effect through the skew scattering or equivalently the 
Mott scattering^iii^ and the side jumpi^iiiii^. The 
skew-scattering contribution diverges in the clean limit 



{(Jxx ~ {e^ /ha){EFT/h) oo) 



oo as 



(6) 



Here, S ^ EsoUimpD/Ep {\S\ <^C 1) is the skewness fac- 
tor with the density of states D, the Fermi energy Ep, 
and the impurity potential strength Wimp, respectively. 
Accordingly, Luttinger reconsidered the issue by means 
of the expansions of axy in Uim^^ and the impurity 
imp'^^, or in h/(EpT) with the relaxation 



concentration n 



imp'^knp^- Then, the leading-order term 



time T ^ h/n 

is proportional to l/(ninipUii„p), which corresponds to 
the skew-scattering contributioit^. The subleading-order 
term, which is of the zeroth order in Uimp and nimp, in- 
cludes the original Karplus-Luttinger's result as well as 
some other terms that partialUy cancel the intrinsic con- 
tribution. The side-jump contribution has the form~ 



(7) 



in the clean limit, with the electron density nd and the 
relativistic Aharonov-Cacher coupling A leading to an en- 
ergy shift by \p X M ■ E due to the applied electric field 
E. It is remarkable that this side-jump contribution is 
insensitive to the relaxation rate, leading to the scaling 
relation p^y oc p^.^. Therefore, it can be incorporated 
into the subleading term in the Luttinger's expansion of 

^xy • 

In the conventional quantum transport theory given 
by Luttinger—, the anomalous Hall conductivity is ex- 
panded in -Eso- Then, the ratio H/{EpT) is the only 
key expansion parameter. The sum of the leading skew- 
scattering contribution and the sub-leading intirnsic and 
other impurity-independent contributions to axy is given 



o2 r 



ha 



+ cEsoD 



(8) 



with c being a constant of the order of unity. Assuming 
= Ep, the ratio EsqD appears only as the overall 



*imp 



factor; 



^xy ~ 7^"^^°^ 



EpT 



(9) 



In this expression, h/ [Epr) is the only relevant parame- 
ter that controls an extrinsic-intrinsic crossover. Namely, 
as far as the expansion in EgoD is valid, the first term 
in Eq. ([H]) as a skew-scattering contribution is dominant 
over the other terms in the clean metal h/{EpT) ^ 1. 
Therefore, it has been believed that the extrinsic skew- 
scattering mechanism is dominan1j2., as supported by 
some experiments around the ferromagnetic Curie tem- 
peratureS and in heavy-fermion compounds showing a 
large susceptibility^i^^. 

Nevertheless, many experimental and theoretical 
works support the Karplus-Luttinger's scenario. At 
a fixed impurity potential strength Uimp, experi- 
mental results on Fe- and Ni-based dilute alloys^, 
CuCr2Se4_a;Br2r^i^, and semiconducting helimagnets 
Fei_yCoj,Si and Fei-j^Mn^SiiS appeared to be consis- 
tent with the Karplus-Luttinger's prediction p^y oc p^^- 
First-principles calculations of the anomalous Hall con- 
ductivity for SrRuOy - and Fe^? in terms of the Karplus- 
Luttinger's scenario or equivalently the TKNN formula^^ 
given by Eq. ([2]) also show large values of <T; 



xy 



which agree with the experimentally observed values at 
low temperatures. The agreement is not consistent with 
a simple-minded perturbation expantion of the intrin- 
sic contribution in ^M, namely, the second term in the 
square bracket of Eq. ([H]) or ([5]) . This points to the impor- 
tance of the resonat enhancement due to the topological 
nature and urges a reexamination of the intrinsic mech- 
anism against the extrinsic one. Theoretically, the two 
characters involved in the anomalous Hall effect have not 
been seriously considered on an equal footing. A unified 
description of both intrinsic and extrinsic contributions 
is called for. It is also helpful to develop a theory of the 
anomalous thermoelectric Hall (Nernst-Eittingshausen) 
effect and the anomalous thermal Hall (Leduc-Righi) ef- 
fect in the same theoretical framework. To provide a 
comprehensive description of the unified theory on these 
anomalous Hall transport phenomena is the main scope 
of the present study. 

This paper is organized as follows: In Sec. |ITJ we 
introduce a model appropriate for studying the inter- 
play between the topological dissipationless Hall cur- 
rent and the extrinsic scattering events. In Sec. IIIIl 
the Keldysh Green's function formalism in the gauge- 
covariant Wigner space is briefiy explained, together with 
the self-consistent T-matrix approximation, the numeri- 
cal results for the equilibrium properties, the Mott rule. 
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and the Wiedemann-Franz law. Then, In Sec. IIVI nu- 
merical results are given for the anomalous electric and 
thermal Hall conductivities, and the scaling relation be- 
tween Gxx and In Sec. |Vl using the Mott rule, we 
show numerical results of the thermoelecrtic Hall con- 
ductivity for the anomalous Nernst effect. In Sec. IVI) 
we clarify relations of the present theory to other theo- 
ries. In Sec. I VIII experimental results on a^y and a^x are 
summarized and compared with the present theory. Con- 
clusions are drawn in Sec. lVIIIl Some necessary details of 
calculations of the Green's functions and the self-energy 
in our formalism are given in Appendices. 



II. MODEL 
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FIG. 2: (Color online) Electronic band dispersions of the 
present model given by Hq. 

A simple model that fully takes into account both the 
"parity anomaly" associated with the avoidcd-crossing 
of band dispersions and the impurity scattering can be 
obtained by expanding the Hamiltonian at a fixed Pz with 
respect to the momentum p measured from the originally 
crossing point po of two dispersions; 



H 
HoiP) 



(10a) 



-Aqo-^ + • (5- X -I- ^(7° (10b) 
2 m 

(10c) 



Will 



30-° ^S{r - rimp) 



with the position r of electron, the Pauli and identity 
matrices & = {d'x,d'y,az) and ct", respectively, and the 
unit vector in the z direction. The first term corre- 
sponds to the level splitting 2Ao = -Eso of two bands at 
the avoided-crossing momentum. The second term gives 
the linear dispersion with the velocity v. The third term 
represents the quadratic dispersion with an effective mass 
m, whose anisotropy has been neglected since it is unim- 
portant. When Aq = 0, Ho{p) given by Eq. (jlObp is 
reduced to the two-dimensional Rashba modeP^ for het- 
rostructures of n-type semiconductors, which possesses a 



magnetic monopole at p = 0. However, note the signifi- 
cant difference in the physical interpretation of the model 
parameters. In the present model for ferromagnetic met- 
als, Aq plays the role of the spin-orbit interaction that 
lifts the accidental degeneracy of two band dispersions 
with the velocity v, while for semiconductors, v is the 
spin-orbit coupling and A is the Zeeman splitting. 

For simplicity, we also assume that the impurity poten- 
tial has the J-functional form of the strength Mimp at ran- 
dom positions {Timp}, since as we will show, this simple 
form bears the extrinsic skew-scattering and side-jump 
contributions. The generalization to long-range potential 
is straitforward but requires rather lengthy calculations. 
Here, note the sgnificant difference in the way of taking 
the impurity average. One might think that the conven- 
tional scheme of taking the average over the impurity po- 
sitions, which we adopt, might be mimicked by taking the 
potential average in the white noise model or the random 
potential model. However, this is true only within the 
Born approximation. In fact, the higher-order scattering 
proccesses, which are crucial for the skew-scattering con- 
tribution, can not be correctly reproduced by the random 
potential model. 

Frist let us consider the model in the absence of the 
impurities. This model has two band dispersions s±{p) 
separated by 2Ap with Ap = ^v^p^ -I- Aq, as shown in 
Fig. [21 Henceforth, the bottom of the lower band in the 
absence of impurities is chosen as the origin of the energy 
and the bottom of the upper band denoted as -Eros = 
£_ {p = 0) is taken as an energy unit. The present model 
possesses the gauge flux 



<7,p 



'2A3 



(11) 



for the band index = ±, as discussed in the literature^. 
Substituting Eq. (fTTj) into Eq. ([2]) and the integration 
over the momentum yields the intrinsic contribution to 
the anomalous Hall conductivity. 



,TKNN 



■'Pit 



when the both bands are partially occupied. Here, p^ is 
the Fermi momentum for the band index a. 

When Ef € [E^s - 2Ao,£;res], t]!^^^ is resonantly 
enhanced and approaches the maximum value e^/2/i. 
Away from this resonance, dominant contributions from 
the momentum region around p = cancel out each 
other or do not appear, leading to a suppression of 
CT™^(« {e^/h){Eso/EF)), and then the perturbation 
expansion in Ego is justified. Therefore, the present 
model, Eq. (I10|) . can be regarded as a minimal contin- 
uum model for a momentum region that gives a major 
contribution to the anomalous Hall effect. 

Actually in the simplest case with the inversion and 
time-reversal symmetry, the massless Dirac fermion 
structure may appear at high symmetry points in pair. 
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They contain Chern numbers with opposite signs as 
in the Honeycomb lattice^'^'i- for graphenes. Then, 
the ferromagnetic moment along the z direction, which 
breaks the time-reversal symmetry, together with the 
spin-orbit interaction avoids the crossing of band disper- 
sions and introduces a energy gap separating two dis- 
persions. This transfers the Chern numbers among dif- 
ferent band indices. In three dimensions, there occurs 
a transfer of the Chern numbers among different values 
of the momentum component p^. In these cases with 
avoided-crossing of band dispersions, a complete inter- 
ference among the topological contribution to the Hall 
conductivity does not occur in general. 

The transfer of the Chern numbers can also occur along 
momentum curves, for instance, along the z direction 
parallel to the magnetization. In this case, the two- 
dimensional massive Dirac-fermion structures studied in 
the present paper appear at each p^, and Aq continu- 
ously changes as a function of Pz with or without a sign 
change. It is likely that there is at least one such struc- 
ture across the Fermi level, which is not accompanied by 
the sign change below the Fermi level. Then, even the 
integration of the two-dimensional anomalous Hall con- 
ductivity over does not lead to a cancellation of the 
topological contribution, and hence it remains of the or- 
der of /ha even for complicated band structures found 
in first-principles calcualtion o^^i^^'^°i'^^'^^ . 



III. QUANTUM TRANSPORT THEORY FOR 
MULTI-BAND SYSTEMS 



We employ the nonequilibrium Green's function 
method based on the Keldysh formalism^^ji^, which has 
recently been reformulated in the gauge-covariant Wigner 
representation for generic multi-component systems^. In 
the linear order in the electromagnetic field, it clarifies a 
systematic way of diagrammatically treating the Smrcka- 
Stfeda formula**^ '"^^ and thus the Kubo formula"*^ with 
the self-energy and the vertex corrections. Imposing the 
self-consistency among the Green's function and the self- 
energy, this automatically satisfies the Ward-Takahashi 
identity. This formalism also reveals that there appear 
two mathematically independent self-consistent equa- 
tions for the linear deviation of quantum distribution 
function in the electric field, one for the Fermi-surface 
contribution and the other for the quantum contribution. 
This is usually not easy to recognize in the integral equa- 
tion for the vertex correction in the Kubo formalism. 



A. Quantum transport theory based on the 
Keldysh formalism in the gauge-covariant Wigner 
space 

Following our previous paper—, we consider the 
Green's functions G" and the self-energies under the 
constant applied electric field which is taken along 



the y direction in this paper, i.e., E = {0,Ey). Here, 
the superscripts a = R, A and < correspond to the 
retarded, the advanced and the lesser components, re- 
spectively. The Green's functions and the self-energies 
are considered in the gauge-covariant Wigner spaced 
composed of the center-of-mass time (T) and space {X) 
coordinates and the mechanical energy e and momen- 
tum p, which can be obtained by the Fourier transform 
of the gauge-covariant derivative. This gauge-covariant 
Wigner representation is advantageous over the other 
choices, since it minimally reduces the arguments of 
the Green's functions and the self-energies. Namely, 
in the uniform steady state, it allows the dependence 
on X^' = {T,X) or X^, = {-T,X) only through the 
electromagnetic potential Af^{X) — {(t){X) , A{X)) or 
Afj_{X) = {-(f>{X), A{X)), which can be totally absorbed 
into the mechanical energy-momentum = (e, p) or 

Then, the Dyson equations are modified with the ap- 
plied field as 

ei- KoiP) ~ti^)\*Gie,p) = 1, (13a) 
Gie,p)-k\el- Koip) -±{e)\ = /. (13b) 

Henceforth, matrices in the Keldysh space are under- 
lined, while those in the band indices are denoted with 
the hat 



G 



So 



G« 


2G< 





G^ 




2S< 







Ho 


M 







a° 







^0 



I 

The symbol 7k- is the Moyal product of the form, 
i{—e)h 



exp 



d pv d pi^ 



(14a) 
(14b) 
(14c) 
(14d) 

(15) 



with the differential operators d and d operating on 
the left-hand and the right-hand sides, respectively, and 
the electromagnetic field tensor F^" = dXf^A'^{X) — 
dx^,A^{X) which is assumed to be constant. The lesser 
Green's function and self-energy play roles of the quan- 
tum distribution function and the vertex correction, re- 
spectively. Since we have assumed that impurity po- 
tential has a (5-functional form, the self-energies are lo- 
cal. The distribution function is now fully quantum- 
mechanical with the e dependence, in sharp contrast to 
the classical Boltzmann transport theory where only its 
integration over e is considered. 

G"(e,p) and S"(e) can be expanded in Ey as 



G"(e,p) 



G^ (e, p) + ehEyG^ {e, p) + 0(i?2), (I6a) 
t'^{e)+ehEyt'^E^{e) + 0{El). (16b) 
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Henceforth, functionals with the subscripts and Ey de- 
note those in the absence of and the gauge-covariant hn- 
ear response to Ey, respectively. Note that even with 
the subscript 0, functionals contain the self-energy orig- 
inating from the impurity scattering. G^ '^ satisfies the 
well-known Dyson equation in the absence of the electric 
field, 



G^-'^{e,p) 



HoiP) 



t^'^{e)]-\ (17) 



' R A <C 

The self-consistent equations for G^' ' are obtained by 

expanding the Moyal product Eq. (fT5|) in the Dyson equa- 
tion in terms of Ey^. It is convenient to decompose 



Gp and Sp into two: 



GXie,p) = GXj{e,p)d,f{e) + GXjj{e,p)f{em 

t<^{e) = t<^j{e)d,f{e) + t<^jj{e)f{e), (19) 

= G^„(£,p)-Gf^(£,p), (20) 

^Xaiie) = ti^{e)-t%{e). (21) 

G^ J and j can be self-consistently determined from 
the quantum Boltzmann equation in the first order in Ey, 



G^ 



^Eyj'^O ^0 ^Ey.I 



^yi Go ~ Gg 



+ ((E^ - t^){d,yG^) + {d,yG^){±^ - ±^)) ,(22) 
or equivalently, 



'^EyJ 



r<R 



id, 



G« 



Py 



Gn 



(23) 



with the velocity Vi(p) — dp.Ho{p), while G^'^ and S^''^ 
are determined from the other self-consistent equation. 



G 



R,A 



^■,R,Af^R.AA,R,A 

^0 ^Ey ^0 



\ (Gf^Vy{d.Gf^) - {dMV)-yGV) . (24) 



Note that the retarded and advanced Green's functions 
are also modified by the electric field, in contrast to 
the single-band case. We also stress that the iJy-linear 
deviation of the lesser component of the self-energy, 
S^^^'*", yields the vertex correction and modifies the p- 
independent current vertex^'*. 

We proceed to the calculation of the self-energies. The 
side-jump and the skew-scattering contributions can ap- 
pear from the first and the second Born amplitudes in 
the vertex correction, respectively^. Furthermore, for 
a large impurity potential strength, impurity states may 
be produced. These nontrivial effects can be fully taken 



into account up to the linear order in the impurity con- 
centration riimp by means of the T-matrix approxima- 
tion. Furthermore, the self-consistent calculation of the 
Green's functions and the self-energies modifies the re- 
sult and allows for producing a nonperturbative effect of 
impurity scattering events, as we will see later. The self- 
consistent T-matrix approximation, which is represented 
by the Feynman diagrams in Fig. [31 gives 



Sg (e) = ninipTg (e) 



(25) 

To'^ie) = Ui„ip ( 1 - Wimp / i^^^y Gf^{e,p) \ (26) 



for the zeroth-order in Ey and 



Sl,7(£) = ni,„pTg^^(e) 



^bI (e) = "imp^g ' (e) 



for the first-order in E^ 



_fp_ 

(27rft)2 

d^p 



G<Ae,p)T^{e) (27) 



G^'^(e,p)rg«'^(£(P8) 



B. Equilibrium properties 

We first solve Eqs. (HZ]), dSS]) and (l26|) for the self- 
consistent T-matrix approximation to obtain the equilib- 
rium Green's functions Gq '^ and self-energies S^'^. All 
the momentum integrations are performed analytically 
for each value of the energy e as given in Appendix IB| 
and then the numerical iteration is repeated until the 
convergence is reached. Figure U] shows thus obtained 
local electron and spin density of states, 



D{s) 
DM 



1 r (Pp 

J (27rft)2 
1 f dPp 



TT J {2nhy 



ImTr 
ImTr 



Gl^{e,p) 



(29) 



G^ie,p)a^ , (30) 



with a set of parameters v = 3.59, Ag 



= 0.1, and 

^/TBoin = 



2mui,„p = 0.2 for three choices of 7Bori 

^^imp^imp 

/h"^ = 0.01, 0.10, and 0.50. Here, 7Bom rep- 
resents the first Born scattering amplitude in the case 
of ti = Ag = 0. The increase of the impurity poten- 
tial strength Uimp lowers the bottom of the band due 
to the broadening of the electron spectral functions and 
also smears out the singularity which is originally present 
at the energy levels of the majority and the minority 
bands at p = 0, i.e., e — e±{p — 0), without the im- 
purities. This broadened spectral feature is beyond the 
semi-classical approximation where the electron spectral 
function has a ^-functional form. This also plays a crucial 
role in eliminating an unphysical singularity in the ther- 
moelectric Hall conductivity a^y at Ep — s±{p — 0) and 
the discontinuity in axy at Ep = e-(p = 0), as shown 
later. 
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FIG. 3; (Color online) Diagrammatic representation of the self-energy in the present self-consistent T-matrix approximation in 
the Keldysh space, which is composed of the infinite series of multiple Born scattering amplitudes. Here, underlined variables 
G and E are matrices in the Keldysh space. Note that the lesser component gives the integral equation for the vertex correction 
in the context of the Kubo formalism. 




butii 



FIG. 4: (Color online) The results of the local electron and 
spin density of states, D(e) and Dz{e), respectively, in the self- 
consistent T-matrix approximation for v — 3.59, Ao = 0.1, 
and 2muimp = 0.2. 



C. Field-induced change of quantum distribution 
function, and the electric conductivity tensor 



Next, the self-consistent results for G^(e,p) and S^(e) 
are plugged into Eqs. (|22|l and (l27l) to calculate G^^ j 
and J self-consistently, and into Eqs. (|24l) and (l28l) 



to obtain the self-consistent solution for G^'"* and E^'"^. 
Details of the calculations of G^^ j and j, and G^'^ 
and S^'^ are given in Appendices [Cl and iDl respectively. 
Then, G| J J is calculated from Ge"^ via Eq. (PU]). 

Finally, the conductivity tensor is obtained as 



^Tot 



<J, 



de 
// 



(2^ 



Tr 



*^{p)G%^{e,p) 



(31) 



with I, j = X, y, where aj°^ has been decomposed into the 
Fermi-surface contribution a^y and the quantum contri- 



7^^ 

'xy 



de f dp 
J {2-Khf 
d^p 



Tr 



ji 



2TTi J {2nh)^^^ 
de f d'^p 



2m J {2TTh)' 



^TTt J {2^h) 



v^{p)G%j{ii,p)\ , (32) 
^iip)G%^ji{e,p)\ fie) 
Vt{p)GE,,iii£,p) 



rTr 



(33) 

The second lines in Eqs. ([221) and ([55)1 are obtained in 
the zero temperature limit. Here, <j(j is completely de- 
termined by the Fermi-surface properties at the zero tem- 
perature, while a(J contains the whole Fermi-sea proper- 
ties and contributes to only the Hall conductivity. Now, 
we can calculate a(j and a(J separately by substitut- 
ing Eq. ^ into Eq. ^ and Eqs. ^ and (EH) into 
Eq. ([55]) . respectively We note that for our practical cal- 
culations, all the momentum integrations in Eqs. (|32p and 
(I33p are performed analytically, and then the remaining 
energy integration is performed numerically in Eq. (|33p . 
Details of the calculations are given in Appendices [Cl and 

El 

Apart from the mathematical separation as given in 
Eq. pip , it is useful to introduce another decomposi- 
tion scheme to distinguish the mechanisms. Effects of 
the scattering events result in the equilibrium self-energy 
Tiq'^ and the vertex corrections associated with E^'"^'^. 
The extrinsic contribution can be ascribed to the ver- 
tex corrections, while the equilibrium self-energy correc- 
tion yields the modification of the equilibrium electronic 
structure like the quasiparticlc dispersion and the damp- 
ing rate. The latter only modifies the intrinsic contribu- 
tion from that obtained for the perfect crystal. Then, it 
is meaningful to separate the total conductivity into the 
intrinsic and extrinsic parts. In particular, the intrinsic 
part is defined as the contribution that survives without 
the vertex correction, 



^Int 



J Int 



ri Int 



(34) 
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where 



J Int 



d£ r. r, f dp 



xTr 



v,{p)G^{e,p)v,{p) {G^{e,p) ~ G^{e,p) ] 
v.{p) (G^{e,p) - G^{e,p)) v,{p)G^{e,p) (35) 



and 



JI Int 



27r- 



dp 



xTT\v,{p)G^{e,p)vy{p)id,G^{e,p)) 

~v,{p){d,G^{e,p))vy{p)G^{e,p) 
~v,(p)G^(e,p)vy{p){d,G^{e,p)) 
+ v,{p){d,G^{e,p))vy{p)G^{e,p)\ (36) 



are obtained from Eqs. (15^ and ([55]) by ignoring 

and E^'"^ in Eqs. (|23p and (l24l) . respectively. The extrin- 
sic contribution is then calculated as the difference 



Ext 
'xy 



Tot 
'xy 



Int 



(37) 



Especially, when the relaxation rate vanishes, a^^ 
analytically expressed asi^ 



/ Int 



IS 



./ Int 



oo I 



-eije- 



dp 



xd,fienip))Imi{np\Vp\n'p) x {n' p\V p\np)) ^ , (38) 

which can be directly derived from Eq. ([22|) or ((23|) . This 
and (T^^ ^"*(t — > oo) compose the Berry-curvature contri- 
bution, i.e., the TKNN formula given by Eq. Namely, 
provided that there is no singular energy dependence in 



the self-energy S^'^, the relation 



TKNN 
'xy 



I Int 

'xy 



oo I 



(39) 



holds, as addressed previouslyi^^. Then, the integration 
by parts shows that ct^^* can be related to the Fermi- 
surface properties^. 

In fact, even with an infinitesimally small impurity 
concentration, i.e., in the clean limit (jiimp 0), cr^^ 
may be suppressed from the TKNN result Eq. ([2|) by 
the contribution from the bottom of the bands due to 
a nearly singular energy dependence of the self-energy 
S^''^, as is shown mathematically in Appendix IB] and is 
plotted in Fig.[4l This reflects that the momentum p and 
thus the Berry curvature are no longer good quantum 
numbers in the presence of the impurity potential. In 
terms of a constant relaxation-rate approximation, a^^y 
can be calculated from the first principle o^^i^^'^°i'^^'^^ . 
Even with a nontrivial energy-dependent self-energy, the 
first-principles calculation is possible, for instance, using 
the GW approximation^. 



Physically, the intrinsic contribution actually corre- 
sponds to axy{i^) in the limit where r — > oo and sub- 
sequently uj Q. In the metallic case, it disagrees with 
the result in the real DC limit (w ^ then t — > oo) 
which is directly relevant to the transport properties. 
When the Fermi level is located within the energy gap, 
crl.y vanishes and al.y agrees with Eq. ((2]). In general, 
al.y is robust against the scattering and thus the vertex 
corrections, and hence we can regard (7^^ as an intrinsic 
contribution even in the presence of impurities. Actu- 
ally, for the present model given by Eq. (|10p . the vertex 
correction to fj^^, namely, the effect of S^. jj on cr^^, 
is canceled out, as shown in Appendix [D] However, the 
Fermi-surface contribution cr^j, is strongly affected by a 
disspation originating from the vertex correction in the 
clean limit, and the extrinsic contribution plays a crucial 
role. 

The expression given by Eq. (|3T|) together with 
Eqs. (|35p and ([36|l coincide with the Smfcka-Stfeda for- 
mulcfc^i^. Moreover, it is remarkable that in the presence 
of scattering, this approach based on Eqs. ((3T|) . ((32)) . and 
(|33l) provides the diagrammatic treatment for the Smfka- 
Stfeda formula^i^, as previously noted^^i^. Here, in- 
stead of diagonalizing the impurity Hamiltonian and ex- 
pressing the conductivity tensor in the diagonalized basis, 
we have taken into account the self-energy and the ver- 
tex corrections due to the impurities, which correspond 
to and , respectively. 

It is also important that formally there exist two inde- 
pendent self-consistent equations for the total quantum 
distribution function (£,p): one is for G^ li^iP) 

and the other for G^''^(e,p), in addition to the 
equilibrium Green's function Gq (£,p) = (G^(£,p) — 
G^(e,p))/(e). Therefore, one can not correctly obtain 
G^ without solving all these self-consistent equations, 
in general. 



D. Mott rule and Wiedemann-Franz law 

Let us consider the electric current J and thermal cur- 
rent J' due to the electric field E and the temperature 
gradient VT applied to the sample. Up to the linear 
order in E and VT, they are conventionally written as 

J^ ^ I](^y^j+«u(-ajT^)), (40a) 

J! = (^"y^^- + ^^^^ (-^j-^)) • (40b) 

i 

Here, Uij is nothing but the electric conducticity ten- 
sor, and aij and the thermoelectric and the ther- 
mal conducticity tensors, respectively. The anomalous 
thermoelectric Hall effect, i.e., the anomalous Nernst- 
Ettingshausen effect, is characterized by a^y 7^ 0, while 
the anomalous thermal Hall effect, i.e., the anomalous 
Luduc-righi effect, by Kxy ^ 0. 
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Smfcka-Stfeda*^ proved that the Mott rule 



3(- 



-T 



(Ty(£) 



6 — ^ 



and the Wiedemann-Franz law 



(41a) 



(41b) 



generaUy hold in the low-temperature limit in the absence 
of inelastic scattering, with the Bolzmann constant fcs, 
the temperature T, and the electric conductivity tensor 
aij (/i) at the chemical potential fi. 

The T-linear coefficient to Hij /T is just proportional to 
Gij, and does not contain new information. On the other 
hand, the low-temperature limit of ctij /T is proportional 
to the derivative of aij {ji) with respect to the chemical 
potential. Therefore, it is calculated from 



Similar techniques have been used to calculate the 
anomalous Hall conductivity in the same model, using 
the Kubo formula and the semi-classical Boltzmann the- 
ory5i. However, whichever method is used, ignoring 
the self-energy correction in the energy integration of 
the equilibrium Green's function in the Kubo formal- 
ism or using the semi-classical distribution function in 
the semi-classical Boltzmann transport theory sometimes 
leads to an unphysical singularity, as we will explain 
in later sections. Unfortunately, this is the case in the 
present model, and one needs to include seriously the 
lifetime broadening of the quasiparticles. In this respect, 
the present fully quantum-mechanical approach gives a 
powerful theoretical formalism by which one can directly 
treat non-trivial quasiparticle spectra modified by the 
self-energy correction, which always eliminates such sin- 
gularity in the presence of a finite scattering strength. 



with 



(42) 



IV. ANOMALOUS ELECTRIC AND THERMAL 
HALL EFFECTS 



a. 



T 



T 



_fp_ 

{2TTh)'' 



■Tr 



(43) 



tTt 



Up)G%ji{ii.p) . (44) 



In this section, we will show the results for the anoma- 
lous electric Hall conductivity axy The low-temperature 
value of the anomalous thermal Hall conductivity K^y can 
be directly obtained from axy with the universal propor- 
tionality constant via the Wiedmann-Frannz law given 
by Eq. (l4Tbl) . 



Note that it includes the Berry-phase curvature at the 
Fermi leve l^^'^^ . 



E. Connection with the semi-classical approach 

Now it is important to clarify the relation between the 
present fully quantum-mechanical approach^^''*'* and the 
semi-classical approaches modified with the Berry phase 
in the momentum spac o^^i^^ . From the present formal- 
ism, if we ignore the self-energy corrections, i.e., E^'"*, 
the transport equations ((22)) and (|24p can be readily in- 
tegrated over the energy e. Note that when we calculate 
the extrinsic transport current from Eq. (j22p . we need 
to multiply the both sides of the equation by the typical 
relaxation rate r to maintain the meaning. Then, they 
are expressed in terms of the semi-classical distribution 
function but of the matrix form. 



Hp) 



de 
27ri 



G^ie,p)-GSie,p))f{e 



Global dependence of extrinsic and intrinsic 
contributions on the Fermi energy and the 
scattering amplitude 



Figure [5] (a) shows the numerical results on the to- 



(45) 



tal anomalous Hall conductivity crj°' 
function of the Fermi energy Ep and the first Born scat- 
tering amplitude 7Born = V'^Born = riimpuf^pm for a 
typical set of parameters, Aq = 0.1, 2muimp = 0.6, and 
2mv^ = 3.59 in an energy unit of i?rcs = 1-0. Henceforth, 
the energy cutoff is taken as = 3.0 and 7Born is varied 
by changing the impurity concentration rtimp which can 
be directly controlled in experiments, while the potential 
strength Uimp is fixed. 

In the clean limit 7Born = f^/TBom 0, crj°* tends to 
increase rapidly in accordance with the extrinsic skew- 
scattering scenario (see Eq. ©). Strength of the diver- 
gence is proportional to Ep in the low electron-density 
limit, and the sign is inverted around Ep = e+(0) = 
ii^rcs ~ SAq. Sign of this skew-scattering contribution 



^skcw 



It is also evident from 



The first term, i.e., the equilibrium part, can be diago- 
nalized by the band representation with the band index 
n and the dispersion e„(p). Then, the energy integration 
yields f{sn{p))- Though this procedure does not diago- 
nalize the second term, i.e., the nonequilibrium part, this 
semi-classical approach usually serves as a good approx- 
imation to the longitudinal conductivity. 



„ also changes by that of Uimp- 
Fig. [5] that a^"" is significantly reduced when the both 
bands are partially occupied, i.e., /x > E^es = e-(p = 0). 
If the self-energy corrections S^'"^ are approximated by 
constant relaxation rates, af^"" might completely van- 
ish when the both bands are partially occupied and a 
discontinuity would also appear when the Fermi level 
crosses the bottom of the minority banctfl. In fact, the 
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FIG. 5: (Color online) (a) The total anomalous Hall con- 
ductivity crjy', (b) the intrisic contribution al^, and (c) the 
extrinsic contribution ct^^' as functions of the Fermi energy 
Ef and the Bron scattering amplitude 7Born in an energy 
unit of Eres = 1-0. The parameters are choisen as v = 3.59, 
Ao = 0.1, and 2m«imp = 0.2. Note the difference of the scales 
for a:cy in (a), (b), and (c). 



self-energy determined self-consistently with the Green's 
function causes a small but finite skew-scattering contri- 
bution even in this case. This is natural since the Fermi 
surface is now not strictly defined and the quasiparticles 
around the Fermi level participate in the extrinsic dissi- 
pative Hall current through the asymmetric scattering. 
We will discuss this issue in detail later in Sec. IVII 

To identify the intrinsic and extrinsic contributions 
quantitatively, we adopt the separation scheme for the 
extrinsic and the intrinsic parts defined via Eqs. (|34[) and 
([57)1 . The intrinsic part a^^^ is plotted in Fig. [3(b) for the 
same set of parameters. Under the resonant condition for 
Ep being around the range \s-{p = 0) — 2Ao, e-{p — 0)], 
a]^y becomes of the order of e^/2/i. With increasing the 
scattering amplitude 7Born, it only gradually decreases 
due to the damping of quasiparticles. Off the resonance, 
CT^"' is significantly reduced by a small factor EsqD. 




FIG. 6: (Color online) The total anomalous Hall conductivity 
(jj"* as a function of Ef and h/r in an energy unit of Ercs = 
1.0. The parameters are choisen as « = 3.59, Ao — 0.1, 
and 2muimp = 0.02. Note the difference of the scale for a^y 
compared with Fig. [5] (a). 



On the other hand, the extrinsic part a^y*' is shown in 
Fig.[5](c). It is evident that the extrinsic skew-scattering 
process always yields a dominant contribution to (Jxy in 
the clean limit. However, with increasing the relaxation 
7Born by the increase of nimp, af^"" rapidly decays and 
becomes comparable to or even smaller than , indicat- 
ing an crossover from the extrinsic skew-scattering regime 
to the intrinsic regime. This rapid decay reflects that the 
skew-scattering contribution originates from intra-band 
processes and hence the skewness factor S remains of the 
order of EsoUimpD/Ep- 

In Fig. [HI we also show the results obtained with the 
same set of parameters except a smaller value of the 
impurity potential strength 2mMinip = 0.02. Compared 
with the case of a larger Uimp shown in Fig. O the skew- 
scattering contribution is smaller, while the intrinsic con- 
tribution almost remains the same. In particular, in the 
level of the (self-consistent) Born approximation instead 
of the (self-consistent) T-matrix approximation, de- 
pends on Uiinp and ni,„p only through 7Born- The reduc- 
tion of the skew-scattering contribution is also natural 
since as explained below Eq. ^ in Sec. [H it is propor- 
tional to l/nimpUiinp (X Uimp/7 for small values of Wimp 
and riimp. Namely, the intrinsic anomalous Hall effect 
becomes more important in this case, and the extrinsic- 
intrinsic crossover becomes clearer. 



B. Extrinsic-intrinsic crossover 

Next, to gain a further insight and to see more clearly 
the extrinsic-intrinsic crossover in the anomalous Hall ef- 
fect, we discuss a^, cr^y'^*^, cr^l and cr^^*, in compar- 
ison with the total value (tJ°* . Figures [7] and [5] show the 
results for several choices of Ep and 2TOUimp as a function 
of the averaged relaxation rate 7 = H/t = 2ImS^^ (Ep) 
at the Fermi level, instead of the Born scattering am- 
plitude 7b. Here, So"(^) is defined by Eq. ((X3|) in Ap- 
pendix 1^ The other parameter values are taken to be 
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FIG. 7: (Color online) ajf, ai'J',criy'"*, ai^, and cr^f at 
i?F = 0.9 as functions of the h/r for the same parameter 
values as Fig. [5] except the impurity potential strength (a) 
2muimp = 0.02, (b) 0.2, (c) 0.4, and (d) 0.6. 



(a) Ep = 0.5, 2muijnp = 0.20 
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(b) Ep = 1.5, 2muinip= 0.20 
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FIG. 8: (Color online) crjf 

a function of the h/r for the same parameter values as Fig. [5] 
with Ef = 0.5, 0.9 and 1.5 for (a), (b), and (c), respectively. 



the same as used in Fig. [Sj 

Let us start with the resonant case with Ep being 
located within the resonant window. We show the re- 
sults for Ep = 0.9 (see the arrow b in Fig. [2]) m Fig. [T] 
It is clear that the extrinsic contribution crj^,^* and thus 
the total Hall conductivity crj°* include the component 
diverging in the clean limit t oo, and the strength 
of the divergence is increased as the impurity potential 
strength 2muimp is varied from (a) 0.02, (b) 0.2, (c) 
0.4, and (d) 0.6. This agrees with the skew-scattering 
scenario a^"" ~ Uia^pEsoDr as given by Eq. As 
noted in Sec. IIV Al the intrinsic contribution as 
well as its components criy'^^ and cr^y are almost un- 
changed by this variation in 2muinip. has a nearly 
saturated value ~ e^/2h when H/t < Aq = 0.1. In- 
creasing the impurity scattering rate "f = H/t from the 
clean limit, cr^y'^^ evolves from 0, reaches the maximum, 
and then gradually decays, while a^y monotonically 
decreases with increasing 7. Accordingly, the total in- 
trinsic contribution is also only gradually decays as a 
function of 7. On the other hand, cr^^ rapidly de- 
cays in proportion to t. Then, as is clear from the 
panels (b)-(d) of Fig. [3 |f7^'Jf| and \a^f\ intersect at a 
value of the relaxation rate h/r proportional to Uimp- 
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This is consistent with the semi-classical arguments: 
'^xy^ ~ ^xy*"" ~ {e'^ / h)ui„ipEsoDT / h becomes compa- 
rable to ~ e^/2/i at h/r ^ Ui^^pEsoD, and even 
smaller than with futher increasing 7 = ^i/r. This 
defines the crossover bewtween the extrinsic anomalous 
Hall effect in the superclean system h/r < UimpEsoD 
and the intrinsic anomalous Hall effect in the moderately 
dirty system fi/r > UimpEsoD. It should also be noticed 
that even within this intrinisc regime, another extrinsic 
contribution is present because of the vertex correction, 
i.e., J, which partially cancels a^Jy and hence reduces 
ctJ°* from , as shown in Fig. [T] 



For a smaU ratio of Ea,o/Ep ^ IQ^^ — 10^^ as in first- 
principles calculations^*'^^ and Ui^pD ^ 1, dominance of 
the intirnsic anomalous Hall effect is realized within the 
usual clean metal of H/t > Ego (several tens of meV). In 
reality, the total Hall conductivity is the sum of the con- 
tributions from all over the Brillouin zone. Since skew- 
scattering contributions from other momentum regions 
are always subject to a similar rapid decay, the above 
extrinsic-to-intrinsic crossover still occurs unless contri- 
butions from all the avoided-crossing regions of band dis- 
persions are mutually canceled out by accident. 
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Next, fixing the potential strength as 2mwimp = 0.2, we 
take different values of the Fermi energy well below and 
above the resonant window. To be explicit, in the panels 
(a) and (b) of Fig. [H we plot the results for Ep = 0.5 
and 1.5 marked with the arrows a and c in Fig. [21 respec- 
tively. In both cases, the intrinsic contribution is much 
reduced from e^/2ft, by a factor of 2Ao/£'rcs- Therefore, 
the perturbation of axy in the spin-orbit interaction en- 
ergy i?so is allowed in these cases. For (a) Ep = 0.5, 
the dominance of the extrinsic skew-scattering contribu- 
tion extends over a wide regime compared with the case 
of Ep = 0.9 shown in Fig. [7] (b). This agrees with ar- 
guments given in the Introduction along the Luttinger's 
theory (see Eq. But with further increase in the 

relaxation rate ^ — h/T, the total Hall conductivity crj"* 
nearly merges into the Fermi-surface part o'^y'^^ of the 
intrinsic contribution. On the contrary to the resonant 
case shown in Fig. [71 not only a^y but also tr^^^^"* is 
finite even in the clean limit. When the Fermi level is lo- 
cated above the resonant window, i.e., Ep > e-{p — 0), 
aly vanishes in the clean limit as shown in Fig. [51 (b) 
for Ep = 1.5, in agreement with the Kubo- formula calcu- 
lation^^, and only gradually evolves into a finite value in 
the presence of a finite damping, which results from the 
self-consistency between the equilibrium Green's function 
G^'^ and self-energy E^''^ in the T-matrix approxima- 
tion, c^j/"' is finite but it suffers from the reduction due 
to the vertex correction associated with j. Then, 

the total Hall conductivity (tJ°' nearly vanishes except 
the skew-scattering contribution. 



FIG. 9: (Color online) Resistivity pxx ~ l/o'xx for the same 
set of parameters as Fig. [5] except a variation of the impurity 
potential strength 2mitinip = 0.02, 0.2, 0.4 and 0.6. (a) and 
(b) are the linear and logarithmic plots. In (b) all the data 
points for 2muimp = 0.02, 0.2, 0.4, and 0.6 are plotted with 
the same symbol for each value oi Ef- 



C. Scaling relations between ctxx and axy 

In the rest of this section, we discuss scaling relations 
between axy and a^x , which is a source of the controversy 
on the interpretation of the experimental results. 

First, we show the results of the resistivity p^x ~ 1 /<^xx 
calculated by using the same parameter values 2mv^ = 
3.59 and Aq = 0.1 in Fig. [1 The data with different 
values of 2mMimp fall into a single curve for each Ep. 
In the clean limit, pxx is proportional to 7 = h/r or 
equivalently <Txx is proportional to r with different co- 
efficients depending on the Fermi energy Ep. However, 
it is clear from both the linear and logarithmic plots in 
(a) and (b), respectively, that for Ep < E^cs, the resis- 
tivity exhibits a different scaling relation pxx oc (/i/r)'^^^ 
with (pxx ~ 0.8 ± 0.05, as in the portion of the Ep — 0.5 
and 0.9 curves. This crossover occurs as the resistivity is 
varied across h/ne'^. 

Figure \TQ\ shows the logarithmic plot of axy against 
axx for the same set of parameters as Fig. [71 In the 
clean limit, the curves nicely follow axy oc axx and the 
ratio axy/axx is proportional to Uimp for a fixed r or 
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FIG. 10: (Color online) Scaling plot of a^y versus (Jxx for the 
same sets of parameter values as in Fig. |8jb) except 2muimp. 
The arrows show the extrnsic-intrinsic crossover scale of a^x 
for 2mMimp ~ 0.20, 0.40, and 0.60. The horizontal solid line 
and the dotted curve represent the values given by the TKNN 
formula of Eq. ((2)| and the intrinsic contribution in the limit 
of 2muiinp — > 0. The vertical dashed line gives the second 
crossover scale of a^xi^ ne^/h) to the 1.6 power-law regime 
in the dirty regime. 



ax'x- Note, however, that for 2muimp = 0.02, the intrin- 
sic anomalous Hall effect is still robust vifithin this range 
of axx, as argued in Sec. IIV Al As axx ^ 2{e'^ /h)EpT/h 
decreases with decrease in r, the relation between axx 
and axy exhibits a upward deviation from the linear one 
axy oc axx , signalling the crossover to the intrinsic regime 
with an almost constant axy As we have already men- 
tioned, this extrinsic-intrinsic crossover occurs around 
o'xx ~ (g^ /ha){Ep /uim-pEa^oD). In terms of the resistiv- 
ity tensor, the crossover occurs when pxx is of the order 
of /xil cm; pxy oc pxx in more conducting region, while 
Pxy oc p^x in Isss conducting region. 

A smaller impurity potential strength Wimp enlarges 
the region of the constant behavior of axy (Note that 
we change riimp to control h/r.) It is remarkable that 
in the case of 2mMimp = 0.02, the intrinsic behavior of 
an almost constant axy is clearly observed in the mod- 
erately diry case. This intrinsic regime with an almost 
constant axy continues even with the variation of axx 
over three orders of magnitude. The magnitude of axy 
in the intrinsic regime is also consistent with experimen- 
tally observed values axy = 10^ — 10'^ f2~^cm~^ in this 
CTajy-constant region, as is summarized in Sec. I VIII 

Further decrease in r changes the scaling behaviour 
to axy oc (axx)^ with ip ~ 1.6. This nontrivial expo- 
nent appears in the dirty regime axx ^ 7re^//i, where the 
longitudinal transport also exhibits the nontrivial scal- 
ing axx oc 7'^'^ and hence the Hall conductivity scales 
as axy oc r'f'^y with ipxy ~ ffxx ^ 1-3. This exponent 
appears irrespective of the position of the Fermi level. 
Actually, the exponent Lp^ has recently been experimen- 
tally verified in various class of ferromagnetic materials, 
as is summarized in Sec. I VIII Note also that this ex- 



ponent approximates to the value expected for the insu- 
lating phase of the quantum Hall system under a strong 
applied magnetic fielci^. These nonperturbative expo- 
nents are obtained as a result of the self-consistency be- 
tween the Green's function and the self-energy and the 
inclusion of the vertex correction. The intuitive under- 
standing requires a further study. It should also be no- 
ticed that the present theory includes neither the weak 
localization corrections steming from the Cooperons and 
diffusons nor effects of the Anderson localization, which 
are crucial to explain recent nontrivial experimental ob- 
servations on disordered thin filma^. 



V. ANOMALOUS NERNST EFFECT 

Before proceeding to discussions and conclusions, we 
give numerical results for the anomalous Nernst effect. 
The T-linear coefficient of the thermoelectric Hall con- 
ductivity Uxy is calculated by evaluating the expressions 
Eqs. (gll), (gll), and (gH) obtained through the Mott rule. 
This includes information on the topological structure of 
the wavefunction at the Fermi level, which can not be 
directly gained from an observation of axy 

First, Ci^°y jT and o;^'^*/'^ for the total and the intrin- 
sic anomalous Nernst effects are shown as a function of 
the Fermi energy Ep and the Born scattering amplitude 
7Born in the panels (a) and (b) of Fig. [TTl respectively. 
Here, the same parameter values as in Fig. [5] have been 
used. In the clean limit, the skew-scattering conrtibu- 
tion o^xy^ yields a dominant and diverging contribution 
for Ep < e+(p = 0) = E^^^ ~ 2Ao = 0.8 (see Fig. [2|), 
reflecting the Ep dependence of a^^y'^ as found in Fig. O 
(a) (or (c)). Besides, there appear two prominent struc- 
tures at Ep = £-|-(p = 0), where the sign change of cr^y™ 
occurs and hence the axy is also strongly enhanced at 
these points with opposite signs. On the other hand, 
oixy is appreciably suppressed in the resonant window 
Ep e [Srcs - 2Ao,£;rcs]- With increasing 7Born, a^.^"'" 
rapidly decays, as in the case of the anomalous Hall ef- 
fect. Figure [TT] (b) shows the intrinsic contribution a]^y 

calculated by imposing the condition of S^'^'^ = 0, as in 

the case of fr^"^*. In the pure case nimp = 0, a]^y is propor- 
tional to the density of the Berry-phase curvature at the 
Fermi level. This shows that the Berry-phase curvature 
is strongly enhanced around Ep — £±{p = 0), and that 
it only gradually decays as a function of the relaxation 
7Born- Then, in moderately dirty systems, there appears 
a crossover from the extrinsic skew-scattering regime to 
the intrinsic regime. However, in the case of the anoma- 
lous Nernst effect, an interference between positive and 
negative contributions appears and a sign change of axy 
often occurs even as a function of 7. Therefore, detailed 
scaling analysis as performed for axy is difficult. 

Note that the divergence of axy/T emerges only in the 
clean limit, in particular, around Ep = £±{p = 0). In 
the clean limit, cr^'y exhibits a kink at Ep = £±(p = 0). 
In the presence of finite relaxation, the quasiparticle spc- 
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FIG. 11: (Color online) (a) The T-linear coefficients of the 
total anomalous thermoelectric Hall conductivity qJ^'/T and 
(b) the intrisic contribution aj."*/r as functions of Ef and H/t 
in an energy unit of E^cs = 1.0. The parameters are choisen 
as u = 3.59, Ao — 0.1, and 2mitinip ~ 0.2, as in Fig. [5] Note 
the difference of the scales for a^y in (a) and (b). 

tra are broadened and this smears out any singularity. 
However, if the broadening effect is ignored, then trj"* 
shows a discontinuity at Ep — £-{p = 0), namely when 
the chemical potential crosses the bottom of the minor- 
ity band, no matter how large is the relaxation associated 
with the impurity scattering^^. According to the Mott 
rule given by Eq. (|41ap . this results in a (5-functional di- 
vergence axy/T — s- oo at this Fermi energy. Therefore, 
one must take into account the finite damping of quasi- 
paticles in performing the energy-momentum integration 
to calculate the conductivity tensor. 

VI. RELATION WITH OTHER THEORIES 

In this section, we clarify the relations of the present 
theory to recent and early theories. 

Recently, the anomalous Hall conductivity in this 
model or similar models have also been intensively dis- 
cussed using the Kubo and/or the Streda formula e^^'^^i^^ 
and the effective semi-classical Boltzmann transport the- 
orj*^. Some of them ignored the vertex corrections to- 
tally^ or partially^'^. Others neglected the self-energy 
corrections or just took a constant relaxation-time ap- 
proximation^i^i^, and further assumed that the relax- 
ation rate is independent of the band index^^, which is 
justified in the limit of Ep £-(p = 0). Then, they 
gave partly different behaviors from our results shown in 



the previous Letter— and the present paper. We stress 
that the vertex corrections are crucial for the transport 
properties and gives a significantly important contribu- 
tion particularly in the clean limit, namely, the skew- 
scattering contribution. In the dilute impurity concentra- 
tion, the scattering potential has a particle-hole asymme- 
try and gives rise to the asymmetric scattering together 
with the spin-orbit interaction for the Bloch electrons, 
even without the potential having a large spin-orbit in- 
teraction studied in the original work^ or /-electronsi^. 
Note that the lack of the self-consistency between the 
equilibrium Green's function and self-energy in the litra- 
tureSiiS is another reason for the discrepancy. 

Now the main disagreement and the controversy on 
the anomalous Hall effect in this model described by 
Eq. (jlOap is on the absence or the existence of the skew- 
scattering contribution and the total anomalous Hall con- 
ductivity in the case where the both bands are partially 
filled, as addressed by Sinova and coworkers^. There 
are two main sources for this discrepancy. In Ref. [5l|, (i) 
the self-energy correction was taken as only two constant 
relaxation rates in the Born approximation, whereas we 
have solved the self-consistent T-matrix approximation 
which becomes exact in the dilute impurity limit, (ii) 
When momentum integrations of the Green's functions 
were performed to calculate axy by means of the Kubo 
formula, they employed the semi-classical approximation 
where the quasipaticle spectrum had a J-functional form 
and hence neglected the self-energy in the denominator. 
Such expansion of axy at the singular point r ^ oo is 
not straightforward to handle in the present model. In 
Ref. Q and the present paper, we have taken the op- 
posite strategy free from the singularity: we start from 
the case with a finite lifetime broadening due to the im- 
purity scattering, and gradually decrease the impurity 
scattering strength. This discrepancy on the anomalous 
Hall effect becomes striking when we consider the anoma- 
lous Nernst effect, as we mentioned in Sec. |Vl Namely, 
applying the Mott rule to the singular results obtained 
for axy in Ref. 'sH, the T-linear coefficient to the ther- 
moelectric Hall conductivity tensor axy diverges when 
the Fermi level crosses the bottom of the upper band, 
i.e., at Ep = S-{p = 0), even with the finite relax- 
ation rate. This is not plausible. From the viewpoints 
of the semi-classical Boltzmann transport theory and the 
Kubo-formula calculation in the Matsubara technique, it 
is required to modify the calculations of these corrections 
beyond the semi-classical approximation. The above ar- 
guments indicate that both the vertex corrections and 
the self-energy corrections are highly important for the 
transport properties and should be properly taken into 
account. 

The present theory explains the anomalous Hall effect 
in the whole regime except in the localized regime. From 
the present results, now the source of the confusion over 
decades is clear. The amplitude of the skew-scattering 
contribution, though it is rather sensitive to details of 
the impurity potential and band structure, can be larger 
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than /h in the superclean case h/r <^ Eso, if we as- 
sume the impurity potential strength of the order of the 
bandwidth or the Fermi energy. In this case, there is no 
chance for the band calculation to reproduce the observed 
value of axy and it is difficult to explain the anomalous 
Hall effect quantitatively from the theoretical viewpoint. 
On the other hand, this skew-scattering contribution de- 
cays for Eso ^5 f^/''' rapidly. The side-jump contribution 
is also small and of the order of {e^ / h){Eso / Ep)^- The- 
fore, the intrinsic contribution, which is of the order of 

/h under the resonant condition, is dominant over a 
wide range of the scattering strength E^o ^ h/r < Ep 
(clean or moderately dirty case). Although Luttinger re- 
considered the Karplus-Luttinger theory^ and gave an 
expansion of axy in Uimp, including the skew-scattering 
contribution as well^ , it fails to reveal the above crossover 
in the space of Ep, Eso Sind h/r. 

Our theory also confirms the condition for the first- 
principles band calculation of the intrinsic anomalous 
Hall conductivity to work reasonably in comparison with 
experiments. It can elucidate the experimentally oberved 
value of axy except a correction arising from a reduc- 
tion due to the vertex correction, when the resonantly 
enhanced intrinsic anomalous Hall effect dominantly de- 
termines the axy This actually occurs in the moderately 
dirty case where axy only weakly depends on the scatter- 
ing rate, as shown in Sec. |TVl 



VII. COMPARISON WITH EXPERIMENTS 

We now turn to the comparison with experimental 
results. The anomalous Hall effect has been investigated 
as a fundamental property in many ferromagnetic ma- 
terials with careful analyses to separate the anomalous 
component from the ordinary one. The results on axx 
and axy at low temperatures are summarized in Fig. [12] 
for Fe, Ni, Co, and Gd films^S, Fe single crystals?", 
SrRuOg^o, Lai_,Sr,CoO3S0, Cui_,Zn,Cr2Se42°, 
Lai_^(SrCa)^Mn03ili^, Nd2(MoNb)207^^ 
Fei_^Co^Sii^, MnSiiS, Gai_^Mn^As^i^>^iS2, 
Ini.^Mn^As^i^, anatase-Co-TiOa^iSiS, and rutile- 
Co-TiO a^^'^" . It is significantly important that all the 
experimental data are categorized into three regimes. 

In the poorly conducting regime, there exists a univer- 
sal scaling relation of axy oc a^x, which agrees fairly 
well with the present theory. This is the case for 
Cui_2;Zn2;Cr2Se4^", Lai-^jSr^^CoOa^, a disordered py- 
rochlore ferromagnet Nd2(Moi_xNba;)207i^, Co-doped 
Mn-doped GaAs^i^. The difference in the am- 
plitudes can be understood as a difference in the num- 
ber of momentum regions with avoided-crossing and/or 
a difference in the relative position of the Fermi level. A 
naive interpretation in terms of axy oc a^^, which can 
be obtained by calculating axy from the perturbative ex- 
pansion in the quasiparticle damping rate, neglecting the 
vertex corrections, and assuming that axx is proportional 
to the damping rate, gives a clearly worse fitting to the 



experimental data than axy oc axx - Though another scal- 
ing behavior of axy oc axx might also be appropriate for 
some experimental data on Gai-^^Mn^jAsi^, axy oc al'^ 
can explain its global dependence. 

In the moderately dirty regime with axx ~ 3 x 10'^- 
5 X 10^ fl^^ cm^^, axy has only a gradual dependence on 
axx and appears to approach constant values of the order 
of 10^-10'^ cm~^ with increasing axx, as observed 
in Fe- and Ni-based dilute alloys^, Cui_a;Zn^Cr2Se4^-, 
SrRuOg^o, metallic foils for Fe, Ni, Co, and Gd^o, and 
MnSi^^. The large amplitude and the robustness against 
the scattering events are consistent with the intrinsic sce- 
nario. The side-jump mechanism also yields an almost 
constant behavior of axy However, it suffers from a 
small factor of Eso/Ep ~ 10~^ — 10~^ coompared with 
the intrinsic Berry-phase contribution under the resonant 
condition, which can not account for the large amplitude 
axy ~ lO'^ cm~^. Therefore, the experimental re- 
sults in this regime can be mostly assigned to the intrinsic 
Berry-phase contribution. Actually, the theoretical curve 
based on the present simplest analysis with Ep = 0.9, 
2mv'^ = 3.59, and 2mwii„p = 0.02, which is shown by the 
red curve, can even explain the gradual dependence of 
axy on axx for Cui_xZn2:Cr2Se4 from the dirty to the 
moderately dirty regime. 

In the superclean regime with axx ^ 5 x 10^ fJ"^ cm~^, 
the curve for axy versus axx tends to deviate from the 
constant behavior, and there appears a rapid increase 
or decrease by another contribution, which is positive or 
negative to the intrinsic contribution, respectively, de- 
pending on meterials. The detailed data seem to depend 
on the properties of dilute impurities embedded in the 
materials. In fact, the ordinary Hall effect is pronounced 
by the Landau-level formation at low magnetic fields 
because of the high mobility. It gives a nonlinear de- 
pendence of the Hall resistivity on the applied magnetic 
field. Furthermore, the remnant magnetization tends to 
decrease with increasing the purity in many highly con- 
ducting materials as Fe and Ni. Therefore, the analysis 
of extrapolating the Hall resistivity to the zero magnetic 
field to obtain the anomalous contribution becomes sub- 
tle, and so far only a few experiments and analyses have 
been performed in this regime. Clearly, futher experi- 
mental studies are required to clarify the scaling behav- 
iors in this regime. 



VIII. CONCLUSIONS 

In conclusions, we have developed a unified theory 
of anomalous Hall effect in ferromagnets, in terms of 
the fully quantum-mechanical transport theory for multi- 
band systems. It confirms that the anomalous Hall effect 
is determined by the intrinsic Berry-phase mechanism 
when (i) the Fermi level is located around an avoided- 
crossing of band dispersions in the momentum space, 
(ii) consequently the magnitude of axy is resonantly en- 
hanced to the order of e^/ (ha) ^ 10'^ ft"-^ cm~^, and (iii) 
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FIG. 12: (Color online) Summary of experimental results on various ferromangets, including transition-metals, perovskite oxides, 
spinels, and magnetic semiconductors. The theoretical curve corresponding to Fig. [10] (2m-!iimp = 0.02) is also shown. The data 
are taken from Miyasato et al^ for Gd, Fe, Ni, and Co films, Fe single crystals, SrRuOa, Lai-^^Sr^^CoOs {x — 0.20, 0.25) and 
Cui-a,Zn^Cr2Se4 {x = 0.0, 0.2, 0.4, 0.5, 0.6, 0.8, 0.9), from Lee et al^ and Lyanda-Geller et al^ for Lai-a.(SrCa)a;Mn03, 
from Iguchi et alM- for Nd2(MoNb)207, from Manyala et alM- for Fei_a;Coa:Si, from Lee et ali^ for MnSi, from Matsukura 
et al^, Edmonds et alM, Yuldashev et al^, and Chiba et alM for Gai-^Mn^^As, from Ohno et alM- and Oiwa et alM for 
Ini_2;Mni;As, from Ueno et al^, Cho et al^, and Ramaneti et al^ for anatase-Co-Ti02 , and from Toyosaki et al^ and 
Higgins et al7° for rutile-Co-Ti02. 



the resistivity p^x is larger than [ha/ e^){E^o / Ep) ^ 1- 
10 /iO cm in Sd transition-metals and 10-100 fj,il cm 
in 5d and rare-earth compounds. In fact, the intrinsic 
contribution suffers from a partial cancellation due to 
the scattering events described as the vertex correction. 
Nevertheless, with these resonant conditions, it remains 
of the order of e'^/h in two dimenions and e^jiha) in 
three dimensions. Then, first-principles calculation can 
give a good estimate of a^y By contrast, in the super- 
clean systems with the lower resistivity, the skew scat- 
tering gives the leading contribution diverging in propor- 
tion to the quasiparticle lifetime. As the damping rate 
increases beyond the energy scale of spin-orbit interac- 
tion, the skew-scattering contribution gradually disap- 
pears and the anomalous Hall effect is dominated by the 
intirnsic contributions which are robust against the relax- 
ation in the conducting regime. This extrinsic-intrinsic 
crossover needs to be verified by further careful exper- 
iments. For dirty systems with axx ^ /ha, another 
new scaling relation <Jxy oc al:^ is obtained. Many ferro- 



magnetic materials are located in this regime, and many 
experimental studies support this scaling relation. The 
present work resolves the long-standing puzzle and con- 
troversy on the mechanism of the anomalous Hall effect 
in a whole region of ferromagnetic metals at low tempera- 
tures, and reveals novel crossovers in quantum transport 
phenomena in multi-band systems. 
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APPENDIX A: NOTATIONS FOR DETAILED 
CALCULATIONS 



First, we introduce the following notations, 

Mp) = 



f-L—0,x,y,z 
f-L—0,x,y,z 



(Al) 
(A2) 
(A3) 

(A4) 



with a = R, A, and <, and 6 = and Ey. Thoughout 
the Appendices, Greek symbols (/i, v, p, ■ ■ ■) are used to 
label the components 0, x, y, and 2, while italic ones (i, 
j, £, • ■ • ) are used only for x, y, and z. Contraction of 
these indices are implicitly assumed in all the expressions 
in the Appendices. 



in the absence of the external fields from the coupled self- 
consistent equations (fT7|) , ([25|) and ([26|) . More explicitly, 
a straightforward calculation show 



I]™(e) = 



^iiiiD^i] 



'T-impMiinpi 



m 



(Bla) 

(Bib) 
(Blc) 

^log^^^44--^-'5o^(s), (Bid) 



0, 



(-Ao + Ef (£)).g„^(£) 
50''" (e) = 



(Ble) 
(Blf) 



APPENDIX B: EXPLICIT FORMS OF G§ AND 

Using the notations introduced in Appendix [Aj it is 
ready to calculate the Green's functions and self-energies with 



is) 



E log (e - pV2™ + h- ^Q°ie) + mv^ + aR^{e)) 



AnhPRi^ie) 

(mw2)2 + 2mv2(e + /i - E™(£)) (-Aq + E«^(e))2 



p=A 



p=0 



(B2) 



(B3) 
(B4) 



For later use, we also introduce Gq through 

G^\e,p) = + S,,{-M + Ef (e))) G^(e,p), 

(B5) 

with the fully anti-symmetric tensor eiji- 

All the momentum integrations are performed analyt- 
ically without any approximation as a function of e, and 
then solve the self-consistent equations for each e. 
APPENDIX C: SOLVING THE 

SELF-CONSISTENT EQUATIONS FOR G^^ j 

AND S| J AND CALCULATION OF ajj 



r 



surface contribution, T,Eyj, self-consistently with that of 
the Green's function, which are defined through Eqs. (fTO|) 
and (|18p . respectively. We exploit the notations defined 
by Eqs. dJT]), and We start with Eq. ((22|) . 

whose component can be explicitly written by using 
ReG^'^ = (G^^-hG^'')/2 and ImG^^ = (G^^-G^^)/2i 
as, 



In the following, we show details in calculating the 
-Ej^-linear deviation of the self-energy due to the Fermi- 
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(Cla) 



ImS^(e)G<°^,(£,p) +ImE^O(£)G<^ + (i?o(p) + ReE^(p)) x G<^j{e,p) 

= ilmG^{e,p))i:<°j{e) + (ImG™(£,p))S<^,,(e) + RcG^(£,p) x ^<^j{s) 
+ z ((IinE™(£))(VpReG^(e,p)) + (ImE^(e))(VpReG™(£,p)) 

- (Vpi?o(p))(ImG™(e,p)) - (Vp7J°(p))(ImG^(e,p)) + (ImE^(£)) x (VpIinG^(e,p))) . (Clb) 
These equations can be expressed in a matrix form 

G<'^^,ie,p) ^ C<^j'\e,p) (KZ{e,p)^<l^,{e)+^Kl^^,{e,p)) (C2) 



with 



,(iJo''(p)+ReEg(e)), (C3) 



(5^JmG^°(e,p) + J,,o(l - (5.o)IniG^''(£,p) + J,o(l - 5^o)ImG^'^(e,p) - e^,pReG^''(£,p), (C4) 
-ImE™(£)(VpReG?"(£,p)) + ImE^(£)(VpReG^(e,p)) 

+ (Vpff°(p))(ImG™(e,p)) - (Vpi?o(p))(ImG^(£,p)) (C5) 
ImE^O(e)(VpReGf - (VpiJ°(p))(ImGf (e,p)) 

+ ImE^Xe)(VpReG«°(£,p)) - (Vpi/^(p))(ImG™(e,p)) 

+ e,,p(ImE^"^(£))(VpIniG^''(e,p)), (C6) 

with the fully antisymmetric tensor e^^p that vanishes if any of /i, and p is 0. Here, the equilibrium properties for 
the self-energy E^'^ with /i = 0, x, y, and z have already been given in Eqs. (jBlap . (jBlbp . and (jBlc^ . and hence those 
for the Green's function G^^ can be obtained from the equilibrium Dyson equation. Then, Eqs. (IC3|) . (|C4[) . (jCSp . 
and (IC6D are rewritten as 



(5pJmE™(£) + (5po^.. + ImSf (e) - e^..(-A/r + ReEf (e)) - v{5f,op, - (C7) 

(5p^ImG^'^(e,p) + (5^0(5^2 + 6f,J^o)lmG^'{e,p) ~ e^^,G§'{e,p), 
-ImE«0(£)(9p^.ReG^°(£,p)) + ImEf (e)(5p^ReGf (£,p)) 
+ fe/m)(ImG™(e,p)) + i;e,,,(ImG^X£,P)), 
ImE™(e)(ap^ReG^Xe,P)) " ip,/m)ilmG^\e,p)) 

+ ImE^Xe)(5p,ReG™(e,p)) + z;e,,,(ImG^"(e,p)) - e,,,(ImEf (£))(9pTmG^^(e,p)), (CIO) 



(C8) 
(C9) 



On the other hand, from Eqs. 



4|, and (jBlfj) . the T- matrix approximation to E^ j, which has been 



(Cll) 



formally described in Eq. (j27[) . is reduced to the following expression in the matrix-form representation. 



''imp "-imp 



2Re ((1 - Uimp5o™(£))"imp5^"(e)) 



2Im ((1 - Uh„pgo™(e))^^imp5o''"(e)) 



(1 - Ui„,pgo™(e)) 



'^fJ.VZ 



.(C12) 



r 



Substituting Eq. (jClip into Eq. (|C2p and integrating over All the momentum integrations are performed analyti- 
p, we obtain cally without any approximation to avoid an unneces- 

d^p 



9e:a') = ^ct::A^) j^c<rie,p)ic^^^^is,p), (cis) 



jT^^o.r('^P)^o%Je,p)B<r{e). 
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sary instability which otherwize occurs in the numerical 
momentum integration, and then the e dependences are 
self-consistently calculated by numerical iteractions. 

Finally, is calculated through Eq. ([5^ . which is 
rewritten as 



Green's functions G'^^(e) given by Eq. p4)) . In the 
following, the integral equations for them are explic- 
itly derived. We also employ the notations defined by 
Eqs. dSIl), (EH), and (1X31) . 



{2nh) 



2 + ^e^fcffl, 



(C15) 



APPENDIX D: SOLVING THE 
SELF-CONSISTENT EQUATIONS FOR Gf^ AND 

Eg AND CALCULATION OF erf/ 

The quantum contribution to the linear response 

77 '^an be calculated from the retarded (advanced) First, each component of Eq. ([M]) is expressed as 



Gi°(£,p) = (Go«°(£,p)^+Go^(£,p)2)Ei°(e) + 2G™(£)Go«(e,p).Eij£) 

- (G^"(£,p)2 - G?(£,p)2)(4E^(£)) • G^{e,p) X (V,i/o(p)) (Dla) 
Gij£,p) = (G„^°(£,p)^ - Go'^le,^)') {^lie) (ds^^is)) x G^(£,p)(VX(p)) 

+ (1 - aeSo«°(£))Go«(e,p) X (Vpi?o(p)) - Go«°(£,p)(aeS«(£)) x (Vp7?o(p)) 

+ 2G?(e,p) (Go«"(£,p)Si°(e) + Go«(e,p)Sij£)) (Dlb) 

Applying the above formula to the present model, we obtain 

Gi°(e,p) = (Go^°(e,p)2 + Go«(e,p)2)l]|0(^) ^ 2G™(e,p)Go«(£,p)Si^(e) - x p{d,J:^^-{s))G§{e,p)' (D2a) 
Gi^^(£,p) = Go«(e,p)Ei^(£) + 2Go«^(£,p)(Gr + Go«(£,p)Sije)) + v^z x p(l - a,I]^«(£))Go«(e,p)2 

(D2b) 

G%{e,p) = Go«(£,p)Eij(e) + 2Gf (e,p)(G„«°(£,p)S|"(e) + G^(£,p)S| (e)) 

+ z; (<5.,((1 - a,I]o^°(£))Go«^(£,p) - (a,I]o«^(e))Gr - (£))Go«(£,p)) Go«(e,p) (D2c) 

Using Eqs. ((26|) and (jBlfj) . components of the self-energy (e) given by Eq. (|28p can be explicitly written as 



'^impufmp ((1 - Wimp.9r(£))' - wLp5?'(e)')~' 

X ai - Wo™(£))' + ^Lp.9o''^(e)') + 2(1 - «i.„pffo™(e))5o''^(£)ffl:(£) 

'^impuLp ((1 " ui^p9o°{e))^ - MLp5?'(e)')~' 

X '{{1 ^i„,p5o™(£))' + ^Lp5o''^(e)') 9§:is) + 2(1 - z.imp5o™(e))ffo''^(£)5f (^) 

^iimp"Lp ((1 - Wimp^rie))" - uLp9o"{£?y\Eli£)- 



(D3a) 

(D3b) 
(D3c) 
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Integrating Eqs. (|D2[) over p and taking into account Eqs. (|B1[) and (|D3p . it is found that 



vS' 



(D4a) 



(27r;i)2 

(fp 
{2nh) 



, Go«(£,p)2 [(1 - 9,I]o™(£))(-Ao + (e)) - (e + ^ - So™(e))(9eSo''^(£))] 



(2^ 



Go«(e,p)2 [(1 - d,J:^%e))i-Ao + (e)) - (e + m - S?"(e))(aeSo''^(£))] 



1 - "impwLp ((1 - Unup9o°ie))^ " uLpg^"(e)2) 



d^p 



{l + vVG^{e,p))G^{e,p) 



with i and being x or y. Then, we obtain 



(2G^°(e,p)I]^;(£) +«(a,So«^(£))Go«(£,p)) G^is,p). 



(D4b) 
(D4c) 

(D5) 



Finally, we obtain 



f{s) [(1 - 9,I]o™(e))(-Ao + (e)) - (e + ^ - sr(e))(a,Sf (e))] 



1 + ?iimpwLp ((1 - Ui,„p5r(e))' - Kr.p9o "{e f) 2 

d^p 



(27r;i)2 

1 - riimpuLp ((1 - Wimpg(f"(£))^ - uLp5^^(e)^) 



(27rS)2 2m 



Gr(£,p)G^'(£,p) 



(l+i;VG^(e,p))G^(e,p) 



(27r;i)2 



-2e-^ftw-^Ini 



de 



(£)) 



d^p p2 
(27r;i)2 2m 



G^{e,pf 



(D6) 



r 



and a trivial result of ct^^ = 0. Note that in the limit of Again, all the momentum integrations are performed an- 

a large momentum cutoff, the square brackets appearing alytically and a special care should be taken in confirming 

in the numerator and the denominator of the first term a convergence of the numerical integration over e for ct^^. 
of Eq. (jD6p coincide, and hence the vertex correction or 



equivalently the effect of S^'"* is canceled out for cr^^ 
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